Patterns of human and bots behaviour on Twitter conversations about sustainability

Sustainability is an issue of worldwide concern. Twitter is one of the most popular social networks, which makes it particularly interesting for exploring opinions and characteristics related to issues of social preoccupation. This paper aims to gain a better understanding of the activity related to sustainability that takes place on twitter. In addition to building a mathematical model to identify account typologies (bot and human users), different behavioural patterns were detected using clustering analysis mainly in the mechanisms of posting tweets and retweets). The model took as explanatory variables, certain characteristics of the user’s profile and her/his activity. A lexicon-based sentiment analysis in the period from 2006 to 2022 was also carried out in conjunction with a keyword study based on centrality metrics. We found that, in both bot and human users, messages showed mostly a positive sentiment. Bots had a higher percentage of neutral messages than human users. With respect to the used keywords certain commonalities but also slight differences between humans and bots were identified.

1. What is the optimal model for identifying account typologies?Can existing models be improved?2. Regarding sustainability, are there differences in the behavioural patterns on Twitter of humans and bots? 3. Are there various patterns of activity?What are their characteristics?4. What is the sentiment on sustainability?Are there differences between humans and bots?How has the sentiment evolved over time? 5. What are the most relevant words used in the tweet text?Are there differences between humans and bots?
In this paper, we present an analysis of the messages posted on Twitter about sustainability with the goal to further understand social opinion on this issue.In particular, the said analysis can accurately identify the account typologies (bots or humans) that at some time made posts related to sustainability.Furthermore, useful information can be collected from knowing interaction patterns and understanding characteristics of these tweets.Specifically, the novelties of this research are: (i) The large database on which the study is built allows us to better understand the patterns of messages published by both humans and bots dealing with a specific topic (in particular sustainability).Notice that previous studies, as the one mentioned above used a much smaller data set.(ii) Building a model based on the digital footprint of each user, which symbolises all types of messages that a user sends, in a similar way as previosly proposed 34 .But in addition to that, we go one step further utilising various compression algorithms (gzip, zlib, bzip2, lzma and smaz) and comparing them.Analogously to 52 , the model is built using attributes related to the user's profile and their activity (which is described through parameters associated with the messages sent by each author (see Sect. "Activity patterns").At the same time, as a novelty, our model utilises these parameters together with others referring to the user's digital footprint.Consequently, a robust model is developed providing optimal results.We apply our model to Twitter data corresponding to users who at some point during the period from 2016 to 2022 sent messages (tweets, retweets, quotes, or replies) related to sustainability.(iii) On the basis of the time series describing the activity of each user, behavioural patterns are detected using clustering analysis and statistical tests.As a novelty with respect to previous research [37][38][39] , the behaviour of humans and bots is studied based on parameters that globally characterise the time series corresponding to tweets, retweets, replies and quotes (see Sect. "Activity patterns").(iv) We determine the most relevant sentiments and words, in addition to detecting some differences according to the activity pattern, the type of account, as well as the polarity of the messages.

Compression algorithms
In this section, based on Twitter data downloaded for the period 2006-2022, we present the main results obtained in our study.The tweets used are those containing the following keywords: sustainable agriculture, sustainable food, renewable energy, green urban, sustainable transport, pollution, sustainable city, and sustainable industry (see Sect. "Building the dataset").We used 38,615 users and 96,252,871 tweets (40,000 users were processed, but only those users with a public account were taken into account in the analysis, which were 38,615).The users for the analysis were randomly selected among all those who at some time sent a tweet with any of the above mentioned keywords related to sustainability.Then, a digital footprint was built for each user (see Sect. "Building the dataset").It consists in a string of characters that represent the different interaction mechanisms in Twitter: post a tweet (' A'), retweet ('T'), reply ('C') and quote ('G') used by the user.As an example, a user's footprint could be a string in the form "ACT CAT TTTAG", meaning that the user has made the following actions: post a tweet, reply, retweet, reply, post a tweet, 4 retweets, post a tweet and quote.
After creating the users' digital footprint a compression algorithm was applied.We applied different compression methods in order to find the most appropiate one to our study.The compression methods used were: (gzip, zlib, bzip2, lzma and smaz).We found that gzip algorithm was the most optimal, see Sect."Building the model").
Figure 1, displays the compression ratio (size of raw digital footprint (in bytes)/size of compressed digital footprint (in bytes)) as a function of the size of the raw footprint for bot and human users (in bytes), both magnitudes are highly relevant 34 for gzip method.It can be seen that this ratio showed a gap between human (dots in blue colour) and bot (dots in red colour) users.For human users with footprints smaller than 3,265 bytes, the compression ratio was in most cases less than 10.However, if the sizes were larger than 3,265 bytes, the ratio showed a large variability.This is in line with the fact that information with high entropy cannot be compressed/ decompressed with optimal efficiency 53 .For bots with footprint sizes lower than 3,265 the compression ratio showed linear growth with size.The linear fit in addition to the raw footprint size that resulted in the highest variability of the compression ratio has been included in the Supplementary Material Document (Tables S1  and S2).It can be observed that bots exhibited a higher median of the raw footprint (see statistical quartiles in Table S3) than human users.

Building the model
We built a supervised learning model based in the use of labelled data sets with the purpose of classifying our used Twitter data to discern whether a user is a bot or not.With the aim of labelling the users as humans or bots, the Botometer tool was utilised (see "Methods" Section).Botometer is a very appropriate resource to identify bots on social networks such as Twitter.However, because bot detection is a complex technique, it has limitations: the accuracy of identification is limited by the quality of the data, and the algorithms used for the analysis.Botometer provides a probability score that a user account is a bot, which can result in false positives and false negatives in detection.As bots become more sophisticated, they more optimally imitate human behaviour.This fact makes it difficult for bot detection tools to keep pace with the changing strategies of malicious users.Human user accounts also exist that tend to be identified as bots.This is because they exhibit repetitive behaviours, or utilise tools to automate certain tasks 54,55 .There is also a linguistic and cultural bias in bot detection algorithms, which implies that tools are not equally effective in all languages 55 .In order to mitigate the limitations of Botometer, we have carried out the following actions: (i) select only tweets in English, (ii) choose only Twitter users who have a public account and (iii) adjust the threshold value for labelling users as bot or non-bot in order to match the proportion of bots values that are indicated in other studies as typical on Twitter.
With respect to the model, it was built using a variety of approaches that correspond to the following three alternatives: • Alternative 1: in this alternative, as input variables, we consider both raw size and compression ratio of the user's fingerprint, similarly to 34 , but additionally we include several compression methods.• Alternative 2: in this alternative, analogously to 52 , certain parameters referring to the messages from an author, as well as to the user profile were considered as input attributes.They were: ratio between amount of persons being followed and number of followers, percentage of retweets, replies, quotes and tweets over total messages.Also included, were the ratio of tweets with multimedia content, and maximum time gap between messages (in hours).No parameters related to compression algorithms were used • Alternative 3: in this option, both alternatives 1 and 2 were taken together.
Various mathematical models were implemented (see Sect. "Building the model"), which were generalised linear (GLM), random forest (RFM) and support vector machine (SVM).They were evaluated utilising several performance metrics.The most optimal mathematical model and compression method were random forest and gzip, respectively.Although, lzma, bz2, zlib also showed a good accuracy (higher than 0.80).For each   S4-S17).The main performance metrics for the three aforementioned alternatives, gzip compression algorithm, and GLM, RFM and SVM procedures are given in Table 1.

Activity patterns
With the purpose to unveil the activity patterns of both humans and bots, a clustering analysis using the K − Means method was implemented for A (tweet), T (retweet), C (reply) and G (quote) interactions.For each user, the following attributes were taken into consideration: mean, standard deviation, median, mode as well as maximum, and minimum number of tweets posted daily, including maximum used lag order and obtained p-value in the Augmented Dickey-Fuller (ADF) test [56][57][58] (see Sect. "Activity patterns").
For the different types of messages analysed, 100 experiments were performed in order to study the interactions of human users.A number of human users identical to the number of bots were randomly selected with replacements in each trial.The clustering tendency was evaluated using the Hopkins statistician (HS).Silhouette and Dunn indexes were also computed to obtain the optimal number of clusters (see Supplementary Material Document, Sect.S2.3).
The results obtained for the different types of interactions are described below.
Type A interactions Data from 3421 bots and 34,229 humans was used in the examination of type A messages.Based on the value provided by both Silhouette and Dunn indexes, the human users could be grouped into two clusters.The average value for HS was 0.01060, with a standard deviation equal to 1.1540 × 10 −05 showing evidence of the existence of clusters (the closer this statistic is to 0, the more evidence in favour of the existence of clusters in the information exist.).The average values corresponding to Silhouette and Dunn indexes were 0.89673 and 0.02408.Taking the 100 experiments into account, the Kruskal-Wallis test 59 showed that there was no difference between groups marked as 1 (group 1), nor between groups marked as 2 (group 2).The average percentage of human users in groups 1 and 2 was 2.24 % and 97.76 % respectively.Regarding bots, according to the Silhouette and Dunn indexes, two clusters were also detected, each having 8.39% and 91.61% of bots.
For human users, Table 2 displays the values of the median of the parameters exhibited by the centroids of each cluster in the 100 experiments performed.Table 2 shows magnitudes corresponding to centroids in each cluster for both humans and bots.

Type C interactions
In order to study the type C messages, 25,264 users (159 bots, and 25,105 humans) were analysed.A variety of the trials executed returned a different optimum number of clusters in accordance with the values of Silhouette and Dunn indexes.Because of this, we were not able to perform an activity pattern analysis.

Type T interactions
In relation to the type T messages, 202 bots and 28,866 humans were explored.The HS, calculated over the 100 experiments, provided an average value equal to 0.06259 with a standard deviation equal to 0.00160.According to the Silhouette and Dunn indexes, 2 was the optimum number of clusters in all trials.On average, each cluster included 7.83% and 92.17% of the total human users.Regarding bots, 3 clusters containing 14%, 4%, and 82 % of the total users were found.
For human users, Table 2, shows the median of the parameters associated with the centroids of each cluster in the 100 experiments.For bot users, the attributes corresponding to the centroids are displayed in each cluster.

Type G interactions
For Type G messages, 3,626 bot and 34,970 human users were examined.Notice that although there is a considerable higher number of bots than in the other types of interactions, the analysis did not detect any cluster either for humans or bots.For human users, the HS showed average and standard deviation values equal to 0.15265 and 0.0011.For bots this statistic was 0.13584.Table 1.Performance metrics of the validation set for the models: generalised linear (GLM), random forest (RFM) and support vector machine (SVM) models, for the three alternatives considered.gzip algorithm has been utilised in Alternative 1 and Alternative 3.

Account typologies
This section aims to investigate whether differences exist between humans and bots behaviour in tweetting.
To do this, we analyze the values of different statistical such as the mean, standard deviation, median, mode, maximum, and minimum number of tweets sent daily, in addition to used lag order and obtained p-value in the ADF test.For A and T messages, the Kruskal-Wallis test demonstrated that considering the above-mentioned factors individually, there were no dissimilarities between the 100 trials performed (p-value > 0.05 ).Therefore, it was possible, without loss of generality, to utilise a single experiment to examine differences with the bot group.For Type A messages, for all considered parameters except the median, a p-value > 0.05 was obtained in the Kruskal-Wallis test proving that a differentiation exists for this magnitude.In relation to Type C, G and T messages a p-value > 0.05 was obtained for all parameters individually analysed.

Sentiment analysis
With the purpose of examining the sentiment of the interactions on sustainability, similarly to 9,11 , only type A messages were taken into consideration, because they are the only ones that bring new personal opinions.
The sentiment from the original tweets posted during the period 2006-2022 was analysed.According to the procedure indicated in Sect."Building the dataset", 40,000 users were selected, but due to the fact that some of them had an empty author field, only 37,650 were considered.For bots and humans, the average text size of each tweet was 114 and 130 characters, respectively.
In the same way to 9,11 for each tweet the sentiment is calculated as an average of the sentiment of the words included in it.The computed sentiment is typified as positive (when its value is in the interval [0, −1]), neutral (if it is 0) and negative (when its value is in the interval [−1, −0]).Only tweets with a subjectivity higher than 0 were considered.
We observed that bot users showed a higher number of neutral tweets than human users.For each user typology and sentiment type (positive, negative and neutral), we built a 17-dimensional vector v tucs , in which each component is the annual average sentiment.For each user typology and period analysed, we create three vectors, describing the positive, negative and neutral sentiments.
Sentiments by user typology, ut, humans or bots: In order to unveil differences between human and bot users behaviour, the cosine similarity between vector pairs corresponding to each sentiment type was calculated.This metric, which presents values in the range [−1, 1], is defined as follows 60 : where a, b symbolise two N-dimensional vectors.a i and b j represent the coordinates of each vector.N is the dimension of each vector.It was observed that, despite the relevant differences in the number of tweets between user typologies (during the analysed period, a maximum 80,000 messages were posted by humans in a year and 14,000 by bots), the cosine similarity was 0.95 for negative and positive sentiment vectors.Neutral sentiment showed a higher difference between both user typologies with a value equal to 0.84. Figure 2 depicts the historical evolution of the sentiment about sustainability (utilising the chosen keywords for the tweets downloaded).
Figure 3 shows the evolution of the average daily sentiment during the first four months of 2022 for both human and bot users.It can be seen that the messages posted about sustainability by bots exhibited a higher range of variability in the daily average sentiment.It varied between a value of 0.001, which is very close to neutrality, up to a value of 0.158.In contrast, the daily average sentiment of human users varied between 0.045 and 0.149.The bots also presented more messages whose average sentiment was located in the aforementioned extreme values (see Supplementary Material Document, Table S18).

Keyword analysis
For each tweet, we obtained the most frequent words, as well as the most relevant ones according to centrality metrics 61 (see Sect. "Keywords analysis").The latter could be assimilated to the process that a human being performs when visualising a text and deduce the most significant content.
For each user typology and cluster, Table 3 shows the top 20 keywords.It can be observed that there are 5 words that are common in all clusters.They are "air pollution" and "new renewable energy" which symbolise important topics in sustainability matters.There are also other common words between clusters, if the results derived from calculating the centrality metrics are compared with those obtained by applying the counting method.It can be seen that some low frequency words are, however, a key element to understand the content of the tweet (see Table 3 and in Supplementary Material Document, Table S20).

Discussion
This research makes a relevant contribution to its application area by extending the work of 34 and 52 , implementing a model that distinguishes between human and bot accounts on Twitter.The incorporation as explanatory variables of the parameters related to footprint compression, profile characteristics and user activity, improves the model's performance.The model was tested with five different compression algorithms (zlib, bzip2, lzma and smaz), showing good results (accuracy > 0.80 ) with four of them slightly favouring gzip as the most optimal.The sample taken into consideration for the construction of the model consisted of 38,615 users with the percentage of bots ranging from 9 to 15 percent 62 .
The time series characterisation associated with each user typology allowed us to perform a clustering analysis for both type A and T messages.For bot and human users, the analysis of type A messages demonstrated the existence of two clusters, pointing out that there were two different ways of operating on Twitter.For bot and human users, type T messages exhibited 2 and 3 clusters, respectively.
The analysis of the evolution of the sentiment of the posted tweets during the period since 2006 to 2022 exhibited interesting patterns.Bots, despite having shown a generally higher percentage of positive posts, exhibited a large number of posts with extreme values, a higher overall polarity and a higher proportion of neutral posts compared to human users.This sentiment characteristic remained relatively stable in the first four months of 2022 for both account types.
If we focus on the keyword analysis, especially that based on the centrality metric, we gain valuable information about the content of the messages.Human clusters show high similarity with each other (0.96).However, the bot clusters, especially cluster 1, exhibited the highest difference with the rest of clusters ( < 0.86 ).In spite of that, five words present in the aforementioned cluster 1 are common to all others.
This research not only expands our understanding of sustainability by analysing patterns of activity, relevant words and sentiment, but also provides practical information for distinguishing between human users and bots on Twitter.Below, we summarise certain limitations of the study carried out in this document, which were previously explained in detail.
Bot detection is a complex task, and the accuracy of bot detection tools is limited by the quality of the data and the algorithms used.In order to mitigate this, in reference to Botometer, several actions were carried out such as selecting only English tweets, choosing only users with public accounts, and establishing a threshold value for labelling in accordance with the values provided by other investigations.With respect to the analysis of sentiments, we use TextBlob due to its simplicity of use, its efficiency, and its processing times that are very appropriate for the volume of data used in this research.
Analysing human and bot activity patterns on Twitter was one of the primary objectives of this investigation, but the implications of the research findings extend beyond this social network.The high importance of sustainability issues in national and international contexts makes understanding the behaviour patterns of human and bot accounts very relevant.This knowledge would enable institutions and authorities, as policy makers and decision makers, to have a more effective influence in achieving more sustainable behaviour.Regarding sustainability matters, companies could also benefit from knowledge of human and bot activity patterns, as well as Table 3. Top keywords related to sustainability got from tweets posted by 37.650 users by cluster and typology (TKG method).Common words to all the clusters are marked in bold.www.nature.com/scientificreports/ the sentiment analysis of messages posted on social networks in order to make their communication strategies more effective.Comprehending how bots can influence opinion on sustainability topics could also be of interest to citizens.All of the above allows us to raise awareness about the need to perform a critical analysis of the information published on social networks.In future research, as a continuation of the investigation described in this document, we plan to implement a comparative sentiment analysis of the tweets utilising different sentiment analysers.Additionally, based on time series, we could implement sentiment prediction models.The procedures to be used could be some including artificial neural networks such as Long Short-Term Memory (LSTM), Adaptive Wavelet Neural Network (AWNN), Elman Recurrent Neural Networks (ERNN) or others such as autoregressive (AR), Autoregressive Integrated Moving Average (ARIMA), in addition to others.Tweets in Spanish could be used to expand the analysis here executed.The influence of bots on the perception of human users could also be explored through contagion models 63 .

Compression algorithms
In order to build a model to be used as one of the explanatory variables, the compression ratio of the raw user fingerprint and several compression algorithms such as gzip 64 , zlib 65 , bzip2 66 , lzma 67 and smaz 68 were used.The gzip method provided the best performance metrics, followed by zlib, bzip2, lzma.smaz algorithm exhibited the worst results, which seems to be because it was specially designed for short texts 69 .It must be noted that a large variability of digital fingerprint sizes are handled in this research.
For the purpose of comparing the methods, gzip, bz2, lzma, smaz, and zlib libraries in PYTHON, were used.They were applied considering the default parameters, in the gzip algorithm the parameters were: mode='rb' (it represents the mode of reading or writing the input/output file, 'r':read, 'b':binary) and compress level= 9 (it can take a value 0 (no compression), or any integer value in the range [1-9].The higher the value, the higher the compression ratio and the lower the processing velocity.The default value is 9.).
zlib 65 method was used considering level = −1 (this attribute can take the value 0 (no compression), -1 (trade-off between velocity and compression ratio), or any integer value in the range [1-9].The higher the value, the more compression, and the slower the computation velocity.The default value is -1 (currently, this value is equivalent to 6)) and wbits = MAXWBITS = 15 (it allow us to manage the window size which was utilised during the compression process as well as whether both header and trailer must have been included in the output.It can take values in the intervals [+9 to +15], [-9 to -15], [+25 to +31].The default value es equal to 15.).
bzip2 66 was applied with mode='rb' and compress level= 9 (these parameters have similar meaning to those explained in gzip algorithm).
lzma 67 was carried taken: mode='rb' (this parameter has an analogous mining to that used in gzip algorithm), check=-1 (it symbolises the integrity type examination to consider in the compressed information.It can take the values: CHECK _ NONE (-1): no integrity check is carried out; CHECK _ CRC32: a 32-bit Cyclic Redundancy Check is performed, CHECK _ CRC64: 64-bit Cyclic Redundancy Check is applied, CHECK _ SHA256: 256-bit Secure Hash Algorithm is carried out.).
smaz 68 algorithm was implemented with mode='r' , and buffering=1 (buffering refers to the mechanism of storing the input data prior to carrying out compression process.This parameter allows us to manage the application of a buffer policy, a value is equal to 1, means that a line is selected as buffering)).
All aforementioned algorithms above are described in the supplementary material document (see Sect.S2.1).

Time series characterisation
For each user, in order to characterise the time series describing the number of daily tweets posted from 2006 to 2022, several parameters were calculated, which were mean, median, mode, standard deviation, maximum and minimum number of messages.Certain data regarding the seasonality of the series were also estimated utilising the Augmented Dickey Fuller test (ADF) 56 .They were maximum used lag order, which was estimated according to 70 , and p-value.It must be noted that with the objective of obtaining correct results in the test, an appropriate selection of lag order had to be carried out 58 .

Building the dataset
In order to build a relevant dataset on sustainability, the keywords to carry out the selection of tweets had to be obtained, to do this task we gathered a small group of persons.If more than 35 keyword candidates were proposed, a clustered and affinity diagram was applied.After this, the obtained keywords were sorted by applying the multiple voting system.Through this procedure, those words that reinforce each other were linked.Finally, the following keywords were chosen for downloading tweets: sustainable agriculture, sustainable food, renewable energy, green urban, sustainable transport, pollution, sustainable city, and sustainable industry.
The reason for reducing the number of words from 35 to 8, as described above, was determined by the keyword selection system itself.The used election method has shown optimal results in other studies 9,11 .The number of words to choose was also conditioned by the feasibility of obtaining, processing and analysing the data.Downloading and processing the tweets from 2006 to 2022 containing these 8 keywords took approximately 12 months, resulting in a volume of 233 GB of raw data and a similar magnitude of processed data.Using 35 words would have lengthened downloading and processing times by involving a much larger volume of data, which could have compromised the novelty of the research.8 is considered an appropriate number of words to achieve an optimal compromise between the amount of information downloaded from Twitter and the processing times needed to conduct the investigation.In the same way to 9 , 10 , and 11 , twarc2 software was utilised to download the Twitter messages.Then, T-Hoarder tool 71 was applied (see Supplementary Material Document, Sect.S1) for formatting the messages provided by Twarc2 from JSON to csv format.The procedure to build the dataset used in this study can be summarized as follow: • 40,000 users were randomly selected among all those who at some time sent a tweet with any of the afore- mentioned keywords.Next, all their messages (tweets, replies, retweets and quotes) from 2006 to 2022 were then downloaded.It must be noted that 2006 was the year in which Twitter was created.From 40,000 users only those with a public account were taken into consideration.In our study we considered 38,615 users that posted 96,252,871 tweets.• In addition to the above, all posted tweets in 2022, including the previously mentioned keywords, were downloaded.684,090 tweets were utilised, which corresponded to 259,170 users.

Preparation of tweets for analysis
For each of the 38,615 selected users, a footprint was built.All messages corresponding to a user were coded according to their type: tweet as ' A' , retweet as 'T' , reply as 'C' and quote as 'G' , resulting in a string of variable size.This footprint can be considered an attribute of each user, similar to human DNA.
In order to build a model that allows us to distinguish between humans and bots, the digital footprint was compressed utilising the algorithms described in Sect."Compression algorithms".
In conjunction with the above, in order to generate an indicator that revealed whether the user was a human or a bot, Botometer API was utilised 72 .This information was used to build the model.Botometer API provided a value in [0,1] range, based on different bot identification methodologies.We assume as bots those users which were characterised with a value > 0.9 .Although other values were tested, this magnitude, as we have already mentioned, provided a percentage of bots in line with that indicated by existing pieces of research 62 .
Additionally, the text of each tweet was also processed.In the same way to 9 , 10 , and 11 various replacements were made: all URLs were changed by the term "LINK", all users' mentions were modified the term "USER", hashtags were modified by term "HASHTAG".Emoticons were transformed to their meaning in text.In addition to the above we also carried out various modifications on the text of the tweet, which were: elimination of conjunctions, terms with a size lower than 3 characters and punctuation symbols.We also changed the figures corresponding to the same term by its lemma.We also corrected the misspelt words and transformed the text to lowercase 9 , 10 , and 11 .The text of the tweets corresponding to 2022, was also identically processed.

Building the model
As explained in Sect."Building the model", in order to detect possible differences between tweets coming from bot and human users, a model to identify the users' categorisation was implemented.The model used several features (according to the selection alternative chosen) and applied various mathematical procedures: Random Forest (RFM), Generalised Linear (GLM) and Support Vector Machine (SVM) algorithms.
Prior to the construction of the model it was necessary to check the correlation between input variables.The Pearson or Spearman method would be applied depending on whether the variables were normally distributed or not, which was examined using the Anderson-Darling test 73 .The following hypotheses were utilised, with a significance level equal to 0.05: (i) H 0 : "The sample derived from a normal distribution" (ii) H a : "The sample did not derive from a normal distribution".
If p-value < 0.05 , H 0 would be rejected.Those variables that exhibited a correlation of more than 0.75 with another variable were eliminated.Similarly to 10,74,75 , in order to build the model, a cross-validation procedure was performed in which NF folds were used.The model was trained NF times, where each time 1 fold was taken as a test set, and each of the remaining NF − 1 folds were utilised as a training set.To estimate the model appropriateness, the mean of an estimated metric (ESTMET) was calculated 75 : ESTMET represents Accuracy, Sensitivity, Specificity, and kappa.A separate end estimation of the aforementioned metrics was also computed utilising the validation sets.NF took a value equal to 5. In order to execute the CV process 80% of the samples were used as training and 20% as validation set.
A short description of the GLM, RFM and SVM has been included in the Supplementary Material Document (Sect.S2.2).

Account typologies
With the aim of describing the activity corresponding to each user, related to tweeting, retweeting, replying, and quoting the following attributes were analysed: • Seasonality of the time series.Parameters related to the Augmented Dickey-Fuller test, which allowed us to examine the existence of this property in the user's activity, these were: (i) Maximum used lag order, www.nature.com/scientificreports/computed as 70 , where several models with different lag orders are fitted and the model with the lowest AIC value is selected.The maximum lag order is defined as the order of the model with the lowest AIC value.(ii) obtained p-value.• For each user, various statistical magnitudes such as number, mean, median, mode, standard deviation, as well as maximum, and minimum number of daily tweets were considered.
In order to examine if there are differences between the above-mentioned variables for humans and bots, the following procedure was performed.
• The normality of distributions of each variable was examined utilising the D' Agostino test 76 .A significance level equal to 0.05 was considered.• If normality existed, the existence of homoscedasticity in each variable should be checked through the Breusch-Pagan test.A significance level equal to 0.05 was utilised.• If normality and homoscedasticity existed, the analysis of variance (ANOVA) method should be executed.
• By contrast, if there was no normality or homoscedasticity in the distributions, the Kruskal-Wallis test would be computed.

Activity patterns
With the goal of detecting behavioural patterns, in both humans and bots, we use time series clustering, which is an unsupervised data mining technique that groups data points based on their similarity.For each type A, C, G, and T messages, the parameters indicated in the previous section were considered.Accordingly, we applied K-Means method to a set of data points SDP: x i , x 2 , x n in which x i = (x i1 , x i2 , ..., x id ) symbolises a vector in R d , and then we get the data distributed into k clusters.Each cluster is represented by a centroid.Given a k value, K-Means method operates as followed 77,11 : • Step 1: In order to represent the first centroids, k data points are picked at random • Step 2: Each data point is assigned to the nearest centroid.
• Step 3: Considering the present cluster memberships, we compute the centroids again • Step 4: In case the convergence criterion is not met, we will repeat steps 2 and 3.
The following convergence criterion is considered 77,11 : • There is no or a very small re-designation of data points to different clusters that exists, or • There is no modification of centroids or there is a minimal change.
• There is a very small reduction in the sum of squared errors (SSE).
• C j Where 77,11 the j th cluster • cen j is the centroid of cluster C j (the mean vector of all the data points in C j ), • d(x, cen j ) is the Euclidean distance between data point x and centroid cen j .
• k is the number of clusters.
For human and bot users, the optimal number of clusters is estimated utilising the Silhouette and Dunn indexes, which are described in the Supplementary Material Document (see Sect.S2.3).

Sentiment analysis
The analysis of sentiment and subjectivity of each tweet was carried out using the Python TextBlob 78 library.This library examines the text of each tweet providing a value for the sentiment and the subjectivity of the tweet in the ranges [-1,1] and [0,1], respectively.This analysis is based on a word lexicon in which each word is described by its positivity /negativity as well as by its subjectivity (subjectivity determines the personal opinion that the text of the tweet contains.The higher the subjectivity, the more personal opinion than objective information there is 9 , 11 ).As we have already explained, the sentiment and subjectivity corresponding to each tweet are calculated by averaging the sentiment and subjectivity for all words included in the tweet.Only English-language tweets were used.
It should be noted that although recent methods such as BERT have proven to be very appropriate for the analysis of complex texts [79][80][81][82][83] , we have used the TextBlob library because in addition to its simplicity of use, it has demonstrated optimal results 9,11,[84][85][86][87] .In this research, as mentioned above, the average sentiment of each tweet was calculated from the sentiment shown for each word in the text of each processed tweet.Due to this and the fact that we handle a large volume of data, with Textblob being faster and less computationally expensive (CPU and GPU) than BERT 88 , we have chosen TextBlob as the sentiment analyser.

Keywords analysis
A keyword analysis was carried out using two procedures, which allowed us to compare the results obtained with each method.www.nature.com/scientificreports/Method 1.In this procedure, once the tweets were processed as indicated in Sect."Preparation of tweets for analysis", a word count was performed, and the most repeated words were considered as the most important.
Method 2. In this approach a similar procedure to the one described in 61 was performed.It is named TKG, which is based on the creation of word graphs, , and examines them utilising certain centrality metrics.
Once the tweets were processed as indicated in Sect."Preparation of tweets for analysis", they were sorted into groups of 30 elements each (this value was taken considering an adequate compromise between network size and computational speed).After this, two types of graphs were constructed from the words included in each tweet.The type 1 graph connected two words, provided one followed the other in the text of tweet TKG 1 .The type 2 graph connected all words that were in the same group of 30 elements ( TKG 2 ).It must be noted that, in the type 1 graph, in order to establish the weights of the links between two nodes i, j, two mechanisms were used: (mechanism 1) the co-occurrence frequency was taken as the link weight, this is, w ij = f ij , and (mechanism 2) the inverse co-occurrence frequency was considered as weight, and this is, w ij = 1/f ij .In the type 2 graph only the mechanism 1 was taken into consideration.
The following centrality metrics were considered: • Degree centrality (DC) of a node i, which is described as: 89 where deg i symbolises the number of links connected to the node i , and • Closeness centrality (CC) of a node i, which is defined as: 90 where V is the set of nodes and d ij is the shortest path between the nodes i and j.

Software programs
Several functionalities were coded in R language, such as: (1) Construction of the model, where the caret and vip packages were used.(2) Cluster and tendency analyses were performed utilising the NbClust and clustertend packages.(3) Execution of statistical tests, in which the nortest and tseries packages were used.(4) Keyword analysis (network building and calculation of centrality metrics), was carried out using the igraph and NLP packages (5) Plotting graphs was implemented using the ggplot2.
Other various funtionalities were also implemented in PYTHON language, which were: (1) Compression algorithms gzip, bzip2, zlib, lzma and smaz, where the packages: gzip, bz2, zlib, lzma and smaz were utilised.(2) Emoji analysis, where emoji and demoji libraries were used to select and replace the emojis by text.(3) Keywords examination, where re package was utilised to correct misspellings as well as to find and extract the hashtags from each tweet.For the removal of irrelevant words, the nltk package was used in conjunction with the stopwords package.For splitting text into tokens and lemmatisation of words, the tokenize and wordnet extensions of the nltk package were used.

Figure 1 .
Figure 1.Compression ratio as a function of the size of the raw footprint for gzip algorithm.Compression ratio is calculated as: size of raw digital footprint (in bytes)/size of compressed digital footprint (in bytes).The arrow points to the critical point of 3265.

Table 2 .
Results of the clustering analysis for the Types of interactions A (post tweets) and T (retweets).For human users, the values of the parameters shown in the table correspond to the median calculated over the 100 experiments.The meaning of the abbreviations is: lag: used lag order in ADF test.H: humans, B: bots, p-value: obtained p-value in the ADF test.mean,sd, median, max and min represent the mean, standard deviation, median, maximum and minimum number of messages.AT account typology, C cluster.When the angle between the two vectors (a and b) is small it generates high cosine values, (close to 1) indicating high cosine similarity.